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Abstract. We present a detailed model describing the effects of wire corrugation on the trapping potential 
experienced by a cloud of atoms above a current carrying micro wire. We calculate the distortion of the 
current distribution due to corrugation and then derive the corresponding roughness in the magnetic field 
above the wire. Scaling laws are derived for the roughness as a function of height above a ribbon shaped 
wire. We also present experimental data on micro wire traps using cold atoms which complement some 
previously published measurements |11| and which demonstrate that wire corrugation can satisfactorily 
explain our observations of atom cloud fragmentation above electroplated gold wires. Finally, we present 
measurements of the corrugation of new wires fabricated by electron beam lithography and evaporation of 
gold. These wires appear to be substantially smoother than electroplated wires. 
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1 Introduction 

Magnetic traps created by current carrying micro wires 
have proven to be a powerful alternative to standard trap- 
ping schemes in experiments with cold atoms and Bose- 
Einstein condensates These so-called "atom chips" 
combine robustness, simplicity and low power consump- 
tion with strong confinement and high flexibility in the 
design of the trapping geometry. Integrated atom optics el- 
ements such as waveguides and atom interferometers have 
been proposed and could possibly be integrated on a sin- 
gle chip using fabrication techniques known from micro- 
electronics. Quantum information processing with a single 
atom in a micro trap has also been proposed 0. 

Real world limitations of atom chip performance are 
thus of great interest. Losses and heating of atoms due 
to thermally exited currents inside conducting materials 
composing the chip were predicted theoretically j8l4j and 
observed experimentally soon after the first experimental 
realizations of atomic micro traps |5I6| . 

An unexpected problem in the use of atom chips was 
the observation of a fragmentation of cold atomic clouds in 
magnetic micro traps |7l8j . Experiments have shown that 
this fragmentation is due to a time independent roughness 
in the magnetic trapping potential created by a distortion 
of the current flow inside the micro wire It has also 
been demonstrated that the amplitude of this roughness 
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increases as the trap center is moved closer to the mi- 
cro wire |1U| . Fragmentation has been observed on atom 
chips built by different micro fabrication processes using 
gold 21] and copper wires |7I8) , and on more macroscopic 
systems based on cylindrical copper wires covered with 
aluminum ^U] and micro machined silver foil ^^l- The 
origin of the current distortion inside the wires causing 
the potential roughness is still not known for every sys- 
tem. 

In a recent letter jll) , we experimentally demonstrated 
that wire edge corrugation explains the observed poten- 
tial roughness (as theoretically proposed in ^21) in at least 
one particular realization of a micro trap. In this paper, we 
will expand on our previous work giving a more detailed 
description of the necessary calculations as well as present- 
ing a more complete set of experimental observations. We 
emphasize that extreme care has to be taken when fab- 
ricating atom chips, and that high quality measurements 
are necessary to evaluate their flatness in the frequency 
range of interest. We will discuss the influence of corru- 
gations both on the edges as well as on the surface of the 
wire and give scaling laws for the important geometrical 
quantities like atom wire separation and wire dimensions. 
We will also present preliminary measurements on wires 
using improved fabrication techniques. 

The paper is organized as follows. In section[21 we give 
a brief introduction to magnetic wire traps and empha- 
size that the potential roughness is created by a spatially 
fluctuating magnetic field component parallel to the wire. 
In section 13 wc give a general framework to calculate the 



2 



T. Schumm et al.: Atom chips in the real world: the effects of wire corrugation 



rough potential created by any current distortion in the 
wire. A detailed calculation of the current flow distortion 
due to edge and surface corrugations on a rectangular wire 
is presented in section 0] In sectional we apply these cal- 
culations to the geometry of a flat wire, widely used in 
experiments. Edge and surface effects are compared for 
different heights above the wire and we present impor- 
tant scaling laws that determine the optimal wire size for 
a given fabrication quality. In section and |7| we show 
measurements of the spectra of edge and surface fluctua- 
tions for two types of wires produced by different micro 
fabrication methods: optical lithography followed by gold 
electroplating and direct electron beam lithography fol- 
lowed by gold evaporation. We also present measurements 
of the rough potential created by a wire of the first type 
using cold trapped Rubidium atoms. 



2 Magnetic micro traps 

The building block of atom chip setups is the so-called 
side wire guide ^ . The magnetic field created by a straight 
current carrying conductor along the z axis combined with 
a homogeneous bias field perpendicular to the wire cre- 
ates a two-dimensional trapping potential along the wire 
(see figure Q). The total magnetic field cancels on a line 
located at a distance x from the wire and atoms in a low 
field seeking state are trapped around this minimum. For 
an infinitely long and thin wire, the trap is located at a dis- 
tance X = fiQ I / {2 TT Bq) . To first order, the magnetic field 
is a linear quadrupole around its minimum. If the atomic 
spin follows adiabatically the direction of the magnetic 
field, the magnetic potential seen by the atoms is pro- 
portional to the magnitude of the magnetic field. Conse- 
quently, the potential of the side wire guide grows linearly 
from zero with a gradient Bq/x as the distance from the 
position of the minimum increases. 

For a straight wire along z, all magnetic field vectors 
are in the {x, y) plane. Three dimensional trapping can be 
obtained by adding a spatially varying magnetic field com- 
ponent Bz along the wire. This can be done by bending the 
wire, so that a magnetic field component along the central 
part of the wire is created using the same current. Alter- 
natively, separate chip wires or even macroscopic coils can 
be used to provide trapping in the third dimension. 

For a realistic description of the potential created by 
a micro wire, its finite size has to be taken into account. 
Because of finite size effects, the magnetic field does not 
diverge but reaches a finite value at the wire surface. For 
a square shaped wire of height and width a carrying a 
current /, the magnetic field saturates at a value propor- 
tional to 1 1 a, the gradient reaches a value proportional to 
I /a?. Assuming a simple model of heat dissipation, where 
one of the wire surfaces is in contact with a heat reservoir 
at constant temperature, one finds the maximal applica- 
ble current to be proportional to a^^^ ^3]. Therefore, the 
maximal gradient that can be achieved is proportional to 
\l\fa. This shows that bringing atoms closer to smaller 
wires carrying smaller currents still increases the magnetic 



confinement, which is the main motivation for miniatur- 
izing the trapping structures. However the magnetic field 
roughness arising from inhomogeneities in the current den- 
sity inside the wire also increases as atoms get closer to 
the wire. This increase of potential roughness may prevent 
the achievement of high confinement since the trap may 
become too corrugated. 

We emphasize that only the z component of the mag- 
netic field is relevant to the potential roughness. A vari- 
ation of the magnetic field in the (x, y) plane will cause 
a negligible displacement of the trap center, whereas a 
varying magnetic field component B^ modifies the longi- 
tudinal trapping potential, creating local minima in the 
overall potential JI] . 



3 Calculation of the rough magnetic field 
created by a distorted current flow in a wire 

In this section, we present a general calculation of the ex- 
tra magnetic field due to distortions in the current flow 
creating the trapping potential. By j we denote the cur- 
rent density that characterizes the distortion in the cur- 
rent flow. The total current density J is equal to the sum 
of j and the undisturbed flow jo e.^ . As the longitudinal 
potential seen by the atoms is proportional to the z com- 
ponent of the magnetic fleld, we restrict our calculation 
to this component. We thus have to determine the x and 
y components of the vector potential A from which the 
magnetic field derives. In the following, we consider the 
Fourier transform of all the quantities of interest along 
the z axis which we define by 

Ai^k{x,y) = / Ai(x,y,z)e^'^'- dz , (1) 

V 2 TT L J 

where we have used the vector potential as an example 
and / stands for x oi y, L being the length of the wire. We 
choose this definition so that the power spectral density of 
a quantity coincides with the mean square of its Fourier 
transform : 

i- j e^^^ {Mz)Am)dz = {\Ai,u?) ■ (2) 

The vector potential satisfies a Poisson equation with 
a source term proportional to the current density in the 
wire. Thus the Fourier component Ai_k satisfies the fol- 
lowing time independent heat equation 

f d d \ 

y-Q^+Q^)^l,k~k'Al^k = -^iojLk■ (3) 

where ji is one component of the current density j. In 
the following, we use cylindrical coordinates defined by 
X = r cos{lp) and y = r shi{ip). Outside the wire, the right 
hand side of equation Q is zero. The solution of this 2D 
heat equation without source term can be expanded in a 
basis of functions with a given "angular momentum" n. 
The radial dependence of the solution is therefore a linear 
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combination of modified Bcsscl functions of the first kind 
/„ and of the second kind Kn- Thus expanding Ai^k on 
this basis, we obtain the foUowing linear combination for 
the vector potential 



(4) 



We retain only the modified Bessel functions of the second 
kind, since the potential has to go to zero as r goes to in- 
finity. The ci^ik) coefficients are imposed by equation (PJ, 
and can be determined using the Green function of the 2D 
heat equation We obtain 



[Ml 



IJkr) 



ji,ki'P,r) rdrdcp . (5) 



Taking the curl of the vector potential and using the re- 
lations K'j-^ — — (/v„_i -I- Kn+i)/2 and 2nKn{u)/u — 
—Kn_i + Kn+i, we obtain the z component of the mag- 
netic field from equation 
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in if 



(6) 

This expression is valid only for r larger than tq, the ra- 
dius of the cylinder that just encloses the wire. At a given 
distance x from the wire, we expect that only fluctuations 
with wavelengths larger or comparable to x contribute to 
the magnetic field, since fluctuations with shorter wave- 
lengths average to zero. Therefore we can simplify expres- 
sion (jnj assuming we calculate the magnetic field above 
the center of the wire {y = 0) for x much larger than rg. 
The argument of In in equation lO is very small in the do- 
main of integration and we can make the approximation 
In{kr) ~ (fcr)"/(2"n!). This shows that the q„ coeffi- 
cients decrease rapidly with n. Keeping only the dominant 
term of the series in equation (O, we obtain 



Bz,k{x) 



^ya 



(fc) 



Ki{kx)] 



(7) 



We will see in the next section that the first factor of this 
expression, characterizing the distortion flow, is propor- 
tional to the power spectral density of the wire corruga- 
tion. The second factor peaks at fc ~ 1.3/x justifying the 
expansion. Fluctuations with a wavelength much smaller 
or much larger than l/x are filtered out and do not con- 
tribute. As we approach the wire, more and more terms 
have to be added in the series of equation 10 to compute 
the magnetic field. We emphasize that the expressions de- 
rived in equations ^ and lO are general for any distorted 
current flow that may arise from bulk inhomogeneities or 
edge and surface corrugations. 




Fig. 1. Rectangular wire considered in this paper. The edge 
roughness and the top surface roughness are illustrated in (c) 
and (b) respectively. 



4 Calculation of the distorted current flow in 
a corrugated wire 

We now turn to the calculation of the distortion in the 
current flow due to wire edge and surface corrugations in 
order to determine the associated c;^ coefficients. We sup- 
pose the wire has a rectangular cross section of width Wq 
and height uq as shown in figure Let us first concen- 
trate on the effect of corrugations of the wire edges, i.e. the 
borders perpendicular to the substrate (model equivalent 
to 23 )■ Figure ^ shows that, in our samples, these fluc- 
tuations arc almost independent of the x coordinate both 
for wires deposited by electrodeposition and by evapora- 
tion. We believe this result to be general for wires fab- 
ricated by a lithographic process, since any defect in the 
mask or in the photoresist is projected all along the height 
of the wire during the fabrication process. Thus, in the 
following, the function f^/i that describes the deviation 
of the right (respectively left) wire edge from ±Wo/2 is 
assumed to depend only on z. 

Conservation of charge and Ohm's law give VJ = 
and J = — X where x is the electrical conductivity and 
V the electrostatic potential. We will make the approxima- 
tion that X is uniform inside the wire. In this case V satis- 
fies the Laplace equation W^V = 0. As we are interested in 
deviations from the mean current density jo = 1/ (uq Wq), 
we introduce the electric potential v = V — joz/x which is 
equal to zero in the absence of deviations. From what we 
have said above, v only depends on y and z and satisfies 
the 2D Laplace equation. The boundary conditions for the 
current density on the wire edge require the current to be 
parallel to the wire edge. Thus v satisfies 



dfr/i 
dz 



(z)x 



dv 

jo~Xg^{y^±Wo/2 + fr/i,z) 
Bv 

-X—{y^±W^/2 + U/uz). (8) 



In the following we assume the amplitude of f^ji to be 
small enough so we can make an expansion to first order 
in /,,/; of both terms. We then obtain a linear relation 
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surface of the wire are 



Fig. 2. Scanning electron microscope images of micro fabri- 
cated wires. Side view (a) and top view (b): electroplated gold 
wire of width 50 fim and height 4.5 fabricated using optical 
lithography. Side view (c): evaporated gold wire of width and 
height 0.7 ^tm fabricated using electron beam lithography. 



between ^(±^^0/2,2) and fr/i{z 
can be written as 



which in Fourier space 



dvk , 



-^—{y = ±Wo/2). (9) 

dy 



The potential v satisfies the 2D Laplace equation, so the 
k component Vk{y) is a linear combination of e^'^^ and 
g-fey fpjjQ coefficients are imposed by the two bound- 
ary conditions of equation To complete the calcula- 
tion of these two coefficients, we introduce the symmetric 
component /+ = (/,, + fi)/2 and antisymmetric compo- 
nent /~ = (/r — /i)/2 of the wire edge fluctuations. Going 
back to the current density, we obtain 



jy,k = ikjo 



cosh(fc y) 
cosh(A:M^o/2) 



sinh(A: y) 
sinh(fc 1^0/2) 



Ik 



(10) 

We note that the symmetric part (first term) of the cur- 
rent deviation is maximal near the wire edges for compo- 
nents with a wave vector large compared to I/Wq. On the 
other hand, the components with a small wave vector are 
constant over the width of the wire. 

We now turn to the calculation of the current distor- 
tions due to surface corrugation. Wc assume the bottom 
surface to be flat, since the wire is supposed to be fabri- 
cated on a flat substrate. We denote by fs the fluctua- 
tions of the height of the wire from its mean value uq (see 
figure We follow the same procedure as for the calcu- 
lation of the effect of the wire edge fluctuations. Now v is 
the electrical potential associated with the current density 
j due to the surface corrugation. It depends on x, y and 
z and satisfies the 3D Laplace equation. To first order in 
fs, the boundary conditions of a current tangent to the 



X^(x = uo, y, z) + jo^— (y, z)^[) 
ox oz 
dv 

= 0,y,z) . 



and 



dv 



±Wo/2, z) = 0. 



(11) 



(12) 



Symmetry arguments show that only the part of fs{y, z) 
which is odd in y contributes to the magnetic field along 
z in the plane y = 0. An even component of fs produces 
currents which are symmetric under inversion with respect 
to the plane y = 0. Therefore, they cannot contribute to 
in this plane. Thus, only the Fourier components 

dy dz 



-ikz 



7rV2L Wo 



sin(2m7r?//Wo) fs{y,z) 



contribute, where m = 1, .., 00. With this definition, fs^. (y) 
2^/tt/Wo J2m=o siii(27rmy/Wo)/sfc „. We choose this defi- 
nition of the Fourier component fs^ „ so that {{fs^ „ P) is 
equal to the 2-dimensional spectral density of fs- To ob- 
tain the electric potential produced by a given component 
fSk m use the expansion 
00 

sin(2m7ry/Wo) = ^ 7™,p sin((2p + Ijny/Wo), (14) 



p=0 



where 



7r. 



-8to 



TT (2(TO-|-p) + l)(2(p-m) + l)' 



(15) 



valid for y G [— VFo/2, VFo/2]. Each p Fourier component 
induces an electrical potential Vk.m,p and, since v satisfies 
the Laplace equation, Vk^m.p is a linear combination of 
P+<^px g^^^ ^-u^x .^^j^j,^^, ^ ^A:2 + ((2p+ l)7r/Wo)2. The 
boundary conditions on the surfaces a; = and x = ito 
determine the coefficients and we obtain 



-XVk,-m,p{x,y) = 
^kjo /5,„7™,p£|^^ Bin((2p + l)7ry/Wo). 



(16) 



With the choice of the expansion (|14|l . the boundary con- 
ditions on y = ztWo/2 arc satisfied by each term. Finally, 
we obtain the current density distribution 



jxk.„, (x, y)=2ikfs^ Joy 

sm}i[vpx) . 



E 

p=0 



Im.p 



sinh(i'pUo) 



sin((2p+I)^y/Wo) 



jvk... (2;, y)=2ikfs^,„Jo Y Wo 

/ cosh(zypa;) {2p + 1)t: 
^ V'"'^sinh(i'pUo) i^pWo 



cos((2p-f l)^2//Wo) 
(17) 



The Fourier components ji^ are obtained by summing 
equation 1)17(1 for m = 1, ... ,00. 
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5 Rough potential of a ribbon shaped wire 

In this section, we combine the results of the two previ- 
ous sections to compute the z component of the rough 
magnetic field in the specific case of a flat rectangular 
wire (ito <SC Wq). This simplification enables us to obtain 
analytical results for a system that is widely used in ex- 
periments l8ll2ll(illVll8liq . 

We do the calculation on the x axis for x > Wq/2 
(and y = 0). Since the wire is considered flat, we replace 
the volume current density j by a surface current density 
a = Jjdx. Then we can rewrite the c/,^ coefficients of 
equation Q as 



Wo II 



dj//n(fcy)[^/.fc(j/)+(-i)"'^i,fc(-y)] 

(18) 

We will first study the effect of wire edge fluctuations and 
give universal behaviors for the magnetic field roughness. 
We will then concentrate on the effect of the top surface 
roughness. We will compare the relative importance of the 
two effects and point out important consequences for the 
design of micro wires. 



5.1 Effect of wire edge roughness 

Let us first study the effect of wire edge fluctuations. Here 
we derive the same results as ^2] in a different way. Note 
that, unlike the calculations in |13) . the calculations pre- 
sented here are only valid for distances from the wire larger 
than (or equal to) Wo/2. The expansion we use is never- 
theless useful because it converges rapidly and permits 
the determination of the magnetic flcld roughness for any 
height larger than VKo/2 after the calculation of a few pa- 
rameters (the c„ coefficients). 

The distorted current flow has no component along 
the X axis. The expression of the rough magnetic fleld is 
then given by the first sum of equation ©. Taking lyS = 0, 
we can rearrange this sum using the equalities Kn{kr) = 
K-n{kr) and Cy_^ = ( — l)"cy„ (see equation ((H^ ). we 
then obtain 

oo 

B,,k = ~k"^{cy,^{k) + Cy,^^,{k))K2n+i{kr). (19) 

n=0 

Since only the Cy^ with even n contribute, we see from 
equation H18|l that only the symmetric part of the current 
density participates to the magnetic field. This is what we 
expect from a simple symmetry argument. For the Cy^^ 
coefficients we obtain 

TT Wq 



L 



^^0/2 , COSh(fcy) 

^^"('^) cosh(fcWo/2) (20) 



As pointed out in the previous section, the sum over the 
angular momenta n in equation l|19(l converges rapidly 



with n if X ^ Wq- More precisely, the dominant term pro- 
portional to Ki{kx) gives the correct result within 10% 
as soon as x > 1.5 Wq- As x approaches x = Wo/2, more 
and more terms contribute, for x = Wo/2, 20 terms have 
to be taken into account to reach the same accuracy. 

We now derive the response function of the magnetic 
field to the wire edge fluctuation for x > Wo/2 which we 
define as R{k,x) — |Sz.A:///^P- '^^ already noticed in 
the previous section, far away from the wire {x 3> Wo), 
only wave vectors k <C 1/Wo are relevant. Then we can 
approximate the integral in equation (|20|l by expanding 
the integrand to zeroth order in k y. Keeping the dominant 
term in the series that defines the magnetic field, we obtain 
the following expression for the response function 



R{k,x) 



47r2a;4 



{kx)^Kf{kx) 



(21) 



For a given height x, as k increases, this function increases 
from zero as fc^, peaks at k = 1.3 /x and finally tends ex- 
ponentially to zero. This behavior can be understood as 
follows. At low wave vectors, the angle between the direc- 
tion of the distorted current flow and the z axis tends to 
zero, thus the contribution of these components becomes 
negligible. At high wave vectors, fluctuations with a wave 
length shorter than the distance to the wire average to 
zero. 

To check the validity of equation (|^ . we plot the di- 
mensionless function i?(fc,a;)/[(/io/)^/(47r^a;'*)] for differ- 
ent ratios x/Wq in figure Q- The limit function corre- 
sponds to a configuration where the distorted current flow 
is concentrated on the line x ^ y = 0. For a smaller dis- 
tance from the wire, the flnite width of the wire becomes 
important and R{k,x) differs from the expression H21|) . 
The amplitude is smaller and the peak is shifted to a lower 
frequency. These effects are due to the fact that as x de- 
creases, the distance to the borders of the wire decreases 
less rapidly than the distance to the central part of the 
wire because of the finite width of the wire. Furthermore, 
because corrugations of high wave vector produce a cur- 
rent density localized near the wire border, their decrease 
in amplitude is more pronounced. 

Assuming a white power spectrum of the wire edge 
corrugations with a spectral density , we can integrate 
the equation pifl over the whole spectral range [20]. We 
then find the following scaling law for the rms fluctuations 
of Bz with the atom-wire distance x: 



X 0.044 



(22) 



This expression is valid for x 3> Wo, the numerical fac- 
tor has been found by a numerical integration of equa- 
tion H21|) . Figure Q shows that this expression is valid 
within 10% as soon as x > 2 Wo. For smaller distances 
X, the fluctuations of magnetic field increase more slowly 
and tends to a constant. The points corresponding to 
X < Wo/2 lie outside the range of the previous calcula- 
tion and their values have been obtained by a numerical 
integration for each x. Note that here J+ is the spectral 
density of /"*". For edges with independent fluctuations. 
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where 




Fig. 3. Response function relating the magnetic field rough- 
ness l-Bzfcl'^ to the wire edge fluctuations (see equa- 
tion 12111 '). Plotted is the dimensionless quantity \Bz^./ f^\'^ x 
47r^2;*/(^oJ)^ as a function oi kx where x is the height above 
the center of the wire (?/ = 0). The different curves correspond 
to different ratios x/Wq going from 0.5 to 4.7 in steps of 0.3. 
Small values of x/Wo correspond to lower curves. The curve 
corresponding to the limit given by equation 1)21^ is also shown 
(dashed line). 
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Fig. 4. Magnetic field fluctuations (Bf) as a function of the 
height above the wire (y = 0). Plotted is the dimensionless 
quantity {B1) Wq /{(fxol)'^ Je), where (B^) is the magnetic field 
roughness and Je is the spectral density of the wire edges as- 
sumed to be white, as a function of x/Wo where x is the height 
above the the wire 1^ . Dashed line: l/i^ law given by equa- 
tion ra . 



J+ = Je/2 where Je is the spectral density of each wire 
edge. The asymptotic behavior of {B^) was first derived 



5.2 Effect of top surface corrugation 

We now consider the effect of corrugations of the top sur- 
face of the wire. As shown in equation H17|l . it induces both 
a current along the x and y direction. The surface current 
densities obtained by integration over x have remarkably 
simple forms. We find 



+ cr. 



(2) 

Vk.n 



'Vk,n 

^(2) 



cos(2^) 



-9ikfc In 27rm / -i y» cosh(fci;) 



sh(fcWo/2) 



(24) 

and K = -^//c^ + (2m7r/Wo)^. In the calculation of a^^ ^ 
the summation over p is not analytical. However, as we 
consider wires with uq <^ Wq, one can make the approxi- 
mation (cosh(i'pMo) — l)/sinh(z/pUo) ~ VpUQ. We then ob- 
tain 



ikfsk.^jo 



—Uq sin(2m7ry/PT/o)- 



(25) 



Comparing equation H25|) and equation (|24|) . we see that 
the current density along x is much smaller than the cur- 
rent density along y provide k, <^ 1/uq (i.e. small wave 
vectors both along y and z). 

Within our flat wire approximation, where only dis- 
tances from the wire x ^ uq arc considered, this is always 
the case. In the following we therefore only consider the 
effect of the current density along y. For x > Wo/2, the 
rough magnetic field is then given by equation p9|l . 

Assuming a white power spectrum for the surface cor- 
rugation of spectral density J5, we now derive some prop- 
erties of the rough magnetic field [201 • For large distances 
above the wire (x ^ Wq), only k components much smaller 
than I/Wq are relevant. Then, as we have already shown, 
the coefficients decrease rapidly with n and the dom- 
inant contribution is given by Cyg . To lowest order in ky, 

Cy„ is proportional to the total current ^^^^i2^v^V^'^y' 

(2) 



J. Thus, the only contribution comes from cTyJ ^ ■ Then, caf 



culations similar to those presented in the previous section 
show that the contribution to {B'l) of the Fourier compo- 
nent m of /sj, is 



{Bi 



J, 



Wq 1 {ii^lf 



X 0.044 



(26) 



where Js is the 2-dimensional spectral density of fs- As 
expected it decreases with m as the contribution of rapidly 
oscillating terms averages to zero for large distances. Com- 
puting the sum over m > gives the scaling law for 
the rms fluctuation of B^ due to surface corrugation with 
atom- wire distance x: 



(Bi 



6 



X 0.044 



(27) 



In figure (0) this expression is compared to numerical cal- 



culations based on equation 1)24(1 . The terms cFyJ^^ con- 
tribute at least 90% of as soon as a; > Wq. 

Comparing edge and surface corrugation, we see that 
for large distances, both effects scale in the same way (see 
equations (|22|l and (|27() ). However, at short distances from 
the wire, the amplitude of the magnetic field roughness 

produced by surface corrugation docs not saturate. In- 

(2) 

(23) deed, although the contribution of crj}^ „ saturates in the 
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effects due to wire edge corrugation, we find: 



Fig. 5. Longitudinal magnetic field fluctuations (B?) pro- 
duced by white noise top surface wire roughness as a function 
of x/Wq (y = 0) )2U| . Plotted is the dimensionless quantity 
{B^)/{Js{iJ.oI)^ /{uqWo)). The dashed lines represents equa- 
tion l|?7jl . 



same way as the effect of wire edge fluctuations, the con- 
tribution of cTy^J^ to the current density diverges as one 
gets closer to the wire. Thus at small distances from the 
wire, we expect surface roughness to become the dominant 
source of magnetic field fluctuations. 



5.3 Consequences for micro wire traps 

The scaling laws (|22|) and H27|l are of major importance as 
they impose strong constrains in the use of micro traps. 
As mentioned in section |21 high magnetic field gradients 
are achieved with small wires and short distances. But as 
the distance to the wire decreases, the roughness in the 
magnetic trapping potential increases. Imposing a maxi- 
mal roughness Z\i?max tolerable in an experiment there- 
fore directly determines the maximal transverse gradient 
accessible with a given realization of a micro wire. 

More precisely, as mentioned in section |21 the maxi- 
mal current in a micro wire is limited by heat dissipation: 

Iinax = S.WoUq ^M- To analyze the scaling of the sys- 
tem, we consider the trap center at a distance comparable 
to the wire width x ~ Wq and a wire height uq small and 
constant. For a given fabrication technology, we expect the 
wire roughness to be independent of the wire dimensions 
Wo and uq and we assume white noise spectral densities 
Je and Js for the edge and top surface corrugations 
Using the above expressions for x and / and equations H22|l 
and (|27|) . we obtain the following scaling laws 



7 2^2 



(28) 



6 Wguo 



for the magnetic field fluctuations induced by the edge 
and the surface roughness respectively. Imposing magnetic 
field fluctuations smaller than Z\_Bniax determines a min- 
imal wire width Wo,min and the maximal transverse gra- 
dient VSmax- If the potential roughness is dominated by 



O.min 



2Tr y 



X 0.044 



1/3 



(29) 



For a potential roughness dominated by effects due to wire 
top surface corrugation, wc find: 



W< 



O.min 



uoAB„ 



27r V Js -5 X 0.044 



1/2 



(30) 



As will be described in the following section, a micro wire 
fabricated by electroplating presents an edge roughness 
of Je — 0.1 //m'^. Assuming a wire without top surface 
roughness, a wire height of uq = 5/im, a typical ^ = 3 x 
lO'^A.m-3/2 and imposing a maximal potential roughness 
of zli^niax ~ IniG, the wirc width is limited to Wo,min — 
700 fim, the maximal gradient will be Z\i?,iiax — 0.2 T/cm. 



6 Probing the rough magnetic potential with 
cold atoms 

In a previous letter [Til, we described measurements of the 
magnetic field roughness produced by a current carrying 
micro fabricated wire. The basic idea is to use the fact that 
the longitudinal density n{z) of atoms along the wire, is 
related to the longitudinal potential seen by the atoms 
through a Boltzniann factor: 



n{z) (X 



,-Viz)/kBT 



(31) 



As discussed in section |2 the potential V{z) is propor- 
tional to the z-component of the magnetic field a the cen- 
ter of the trapping potential. Our typical thermal energy, 
1 /iK, corresponds to a magnetic field of 15 mG for a ^^Rb 
atom in the F = 2,mF = 2 state. Since longitudinal den- 
sity variations of order 10% are easily visible in our exper- 
iment, we are sensitive to variations in the magnetic field 
at the mG level. 

The micro wire we used to create the magnetic po- 
tential is a 50 /im wide electroplated gold wire of 4.5 ^m 
height (see figure (0)). The process of micro fabrication is 
the following: a silicon wafer is first covered by a 200 nm 
silicon dioxide layer using thermal oxidation. Next, seed 
layers of titanium (20 nm) and gold (200 nm) are evap- 
orated. The wire pattern is imprinted on a 6 fim thick 
photoresist using optical UV lithography. Gold is electro- 
plated between the resist walls using the first gold layer 
as an electrode. The photoresist is then removed, as well 
as the gold and titanium seed layers. Finally the wire is 
covered with a 10 /im layer of BOB resin and a 200 nm 
thick layer of evaporated gold. The gold surface acts as a 
mirror for a magneto-optical trap. The procedure for de- 
ducing the potential roughness from images of the atomic 
cloud is complex and wc refer the reader to 
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Fig. 6. Rough magnetic field Bz{z) normalized to the cur- 
rent in the micro wire. Solid lines: magnetic field measured 
using cold atomic clouds. Dashed lines: magnetic field calcu- 
lated from the measured corrugation of the edges of the wire. 
The different curves have been vertically shifted by 0.1 G/ A 
from each other and heights above the wire are indicated on 
the right. 



Figure lO shows the measured longitudinal potential 
for various distances above the wire. We also show the 
power spectral density of these potentials in figure O . A 
region of 1.6 mm along the wire is explored by the atoms. 
To estimate the power spectral density of the potential 
roughness we divide the total window in three smaller 
windows overlapping by 50% \21\ . In each window, the 
fourier transform of the potential is computed after mul- 
tiplication with a Hamming window and the estimate of 
the spectral density is the average of the square of the 
Fourier transforms. 

In figure 10 , a flat plateau is visible at the highest wave 
vectors {e.g. k > 0.07 /xm~^ at 46 /xm and k > 0.04 /im~^ 
at 80 fim). The level of this plateau depends on experi- 
mental parameters such as the temperature and density of 
the atom cloud. On the other hand the spectral density at 
low wave vectors, i.e. in the region where it rises above the 
plateau, is independent of these parameters. This obser- 
vation leads us to conclude that while the low wave vector 
part of the spectrum corresponds to a potential seen by 
the atoms, the plateau at high wave vectors is due to in- 
strumental noise in our imaging system, such as fringes. 
We expect it to vary in a complex way with tempera- 
ture and atom density. Qualitatively, smaller atom-wire 
distances, which are analyzed with higher temperature 
clouds, should result in higher plateaus. This tendency 
is indeed observed in figure IJJ). 

To measure the wire corrugations, we removed the 
atom chip from vacuum and etched off the gold mirror 
and the BCB layer. We analyzed the bare wire with scan- 
ning electron microscopy (SEM) and with atomic force 
microscopy (AFM) techniques. The function / describ- 
ing the edge corrugation is extracted from SEM images 
such as ^p). Rms deviations of the edges are as small 
as 200 nm, and we use a 50 fiia x 50 fim field of view 
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Fig. 7. Spectral density of the magnetic field roughness for 
different heights above the wire. The points represent experi- 
mental data. The curves result from the calculations detailed in 
the text. Solid curves: expected noise due to wire edge rough- 
ness. We used the power law fit to the spectral density of the 
wire border fluctuations. Dashed curves: expected noise due to 
top surface roughness. 



in order to have a sufficient resolution. We use 66 over- 
lapping images to reconstruct both wire edges over the 
whole wire length of 2.8 mm. We identify no correlation 
between the two edges. The spectral density obtained for 
/+ = (/i + /r)/2 is plotted in figure ©. We see two 
structures in the spectrum: first, we observe fluctuations 
with a correlation length of 0.2 /im and 100 nm rms ampli- 
tude. It corresponds to the fluctuations seen on figure ijSb) 
which are probably due to the electrodeposition process. 
Second, roughness with low wave vectors is present and 
raises significantly the power spectral density in the 0.01- 
0.1 /im~^ range. For the spectral range 0.01-1 ^m~^, the 
wire border fluctuations are well fitted by a power law 
J = 3.2 X lO-^fc-2-15 4- 8.2 X W^^fim^ as seen in fig- 
ure (jSJ. We use this expression to compute the spectra 
shown in figure (0. As we measured / over the whole re- 
gion explored by the atoms, we can not only compare the 
spectral densities of the magnetic field roughness but we 
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also can compare the direct shape of the magnetic field 
Bz{z). This is done in figure @ where the magnetic field, 
computed from / as described in the previous sections, is 
shown by dashed lines. We note that no adjustment has 
been applied to superimpose the two curves, the absolute 
position of the atoms with respect to the wire is known to 
the 3 /im resolution of our imaging system. 

A different approach is possible to check consistency 
between the wire edge measurements and the potential 
roughness measured with cold atoms. As seen in equa- 
tions idni) and 12011, |-Bz,fc(a;)P = J+i?(fc,a;), where the re- 
sponse function R{k, x) does not depend on the wire edge 
roughness. We compute R{k,x) using the expansion l(T^ 
and for a given k component, we deduce J+ by fitting 
the decay of \Bz,k{x)\'^ with height (see figure O). In fig- 
ure JHJ, the values of J+ obtained by such a procedure are 
compared to the function measured with the electron 
microscope. We find good agreement. 

The corrugation of the top surface of the wire is mea- 
sured using an AFM and the observed power spectral 
density is plotted in figure The spectrum is flat for 
wave vectors smaller than 1 ^m~^ with a value Js = 
1.6 X 10^^ yum*. Unfortimately, we were not able to obtain 
the spectrum for very long wave vectors. For purposes of 
calculation, we shall simply assume that the spectral den- 
sity below 0.1 fim~^ has the same value as between 0.1 and 
1 ^m~^. The result of this calculation is plotted figure iQ 
(dashed lines). 

Our results indicate that the magnetic field roughness 
measured with cold atoms is explained by wire corruga- 
tion. At low wave vectors (fc < 0.04 ^m""'^), it seems that 
the magnetic field roughness is primarily due to edge cor- 
rugations. The good agreement between the observed field 
and the calculation shown in figure l|njl are the strongest 
evidence for this conclusion. For wave vectors larger than 
about 0.5/im~^, the corrugations of the top surface are 
expected to contribute as strongly as those of the edges. 
This wave vector regime however, is not being stringently 
tested by our data. Since we have no data on surface corru- 
gation at wave vectors below 0.1 ^m^^, it is possible that 
the contribution from this effect is larger than shown in 
figure 0. The atom data in the figure however, indicate 
that the surface effect is not the dominant one although 
given our signal-to-noise it could be of comparable mag- 
nitude. 

Figure also suggests that there might be a grain 
structure in the bulk of the wires: we have no additional in- 
formation on possible current deviations due to this effect, 
but the success of the model based on wire edge roughness 
seems to indicate, that it is not important in our system. 



7 Improved fabrication process for micro 
wires 

The fabrication technology described above limits us to 
atom wire separations greater than several tens of mi- 
crons if we want to obtain a reasonably smooth potential. 
In order to improve the quality of our wire, we turn to a 
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Fig. 8. Measured spectral density of the edge roughness of the 
electroplated wire (upper curve) and of the evaporated wire 
(lower curve). For the electroplated wire, the spectral density 
of /"'" = (/i -|- /r)/2 is plotted. For the electroplated wire, Jf /2 
is plotted, where J/ is the spectral density of a single border 
of the wire, as expected for the spectral density of for in- 
dependent wire border fluctuations. The crosses indicate edge 
roughness for different k vectors reconstructed from the decay 
of the corresponding roughness measured with the atoms. The 
solid line is a polynomial fit to characterize the edge roughness 
in the region of interest. 




Fig. 9. Power spectral density of the wire top surface rough- 
ness measured with an AFM. We plot the spectral density cor- 
responding to the transverse mode m = 1 (ky = 2Tr/Wo) which 
is the first one to contribute to magnetic field roughness. The 
horizontal line indicates the mean value for fcz ranging from 
0.2/im~^ to lfim~^. 



different micro fabrication process similar to |22| : the wire 
structures are patterned onto an oxidized silicon wafer us- 
ing electron beam lithography. We use gold evaporation 
and a standard lift-off technique to obtain 700 nm square 
cross section wires as shown in figure (|2;). 

We extract the wire border roughness from SEM im- 
ages and the obtained power spectral density is plotted in 
figure 10) (lower curve). In the spectral range studied, the 
roughness is greatly reduced compared to the first fabrica- 
tion process. This was expected as the grain size of evapo- 
rated gold is much smaller than of electroplated gold. Un- 
fortunately, we do not have a quantitative measurement of 
the power spectral density in the 0.01-0.1 /im^^ range. In- 
deed, as we had to reduce the field of view to increase the 
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resolution, it becomes very difficult to overlap hundreds 
of SEM pictures without adding spectral components due 
to stitching errors. We still hope to also have reduced the 
wire edge roughness in this frequency domain, as it has 
been demonstrated recently by the Heidelberg group [23] 
using similar wires. 

Gold evaporation produces surfaces of optical quality 
at visible light. Thus the roughness of the top surface of 
the evaporated wire is expected to be much smaller than 
that of an electroplated wire. 



8 Conclusion 

Our goal in this paper has been to give a more detailed de- 
scription of the work which led to our conclusion that wire 
corrugations can account for the magnetic field roughness 
typically observed in atom chip experiments. We wish to 
emphasize in this paper that great care must be taken to 
characterize the roughness of a micro fabricated wire. The 
ratio of the rms roughness to the wavelength of the im- 
perfections is below 10^^. Thus a single microscope image 
cannot reveal the imperfections. 

The model we use has already been suggested in ref- 
erence • Here we have given more details of the calcu- 
lation as well as some physical arguments explaining the 
results. We have also extended the calculation to include 
the effects of corrugations of the top surface of the wire. 
The top surface corrugations become increasingly impor- 
tant as the distance to the wire decreases, while the effect 
due to wire edge roughness saturates. 

The equations (|^ and lf?7jl , giving the behavior of the 
magnetic field roughness due to edge and surface corru- 
gation as a function of height, are important scaling laws 
that one should keep in mind in the design of atom chips. 
The requirements of small roughness and high transverse 
confinement impose a tradeoff in choosing a wire size for 
a given fabrication quality. We do not believe however 
that we are at the end of our progress in improving the 
fabrication technology. Thus sub-micron scale atom chips 
continue to hold out much promise for the manipulation 
of ultra cold atoms. 

We thank David Hermann for help in calculations. 
This work was supported by The E.U. under grant (IST- 
2001-38863 and MRTN-CT-2003-505032), as well as by 
the DGA (03.34.033). 
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